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ABSTRACT 

A  comprehensive  analysis  of  the  available  literature  on 
the  atmosphere  of  Jupiter  is  made,  and  two  models  are  selected 
as  being  most  likely.  It  is  suggested  that  a  simple  probe  to 
Jupiter  might  resolve  the  question  of  which  atmosphere  actually 
exists . 

The  general  ballistic  entry  equations  are  displayed,  and 
specific  solutions  for  direct  and  displaced  orbital  entries 
are  made  in  two  dimensions  for  Jupiter.  It  is  found  that  a 
large  dense  vehicle,  requiring  the  Saturn  booster  and  utilizing 
a  graphite  ablative  heat  shield,  can  make  a  direct  parabolic 
entry  at  angles  up  to  30°  •  It  is  doubtful  that  mid-course  and 
terminal  guidance  exists  which  can  give  the  precise  injection 
angles  required  by  this  mission. 

The  feasibility  of  atmospheric  determination  by  measure- 
ments of  drag  load  factor,  altitude,  and  heating  rate  is 
discussed.  The  general  conclusion  is  that  if  a  successful 
entry  can  be  made,  the  atmosphere  can  be  determined.  The 
success  of  the  entry  depends  on  control  of  the  entry  angle. 
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OBJECT 

The  object  of  this  thesis  is  to  examine  the  vehicle  design 
parameters  for  an  entry  into  the  atmosphere  of  Jupiter,  and  to 
investigate  the  feasibility  of  atmospheric  determination  by  a 
simple  entry  vehicle. 


CHAPTER  1 
THE  ATMOSPHERE  OF  JUPITER 


1.1  Introduction 

The  purpose  of  this  thesis  is  to  examine  the  design  parameters 
for  a  ballistic  entry  vehicle  into  the  atmosphere  of  Jupiter.  The 
first  step  is  to  determine  the  most  probable  model  of  the  Jovian 
atmosphere. 

The  Jovian  atmospheres  are  quite  different  from  those  of  the 
Terrestrial  planets.  Jupiter  is  the  closest  and  easiest  to  reach 
of  these  major  planets,  and  therefore  will  provide  Man's  first 
experiment  with  such  an  atmosphere. 

The  most  generally  accepted  theory  of  planetary  formation  is 
Kuiper's  protoplanet  theory  (ref  l).  The  Solar  nebula  was  broken 
up  by  gravitational  instability  and  turbulence  into  the  several 
protoplanets  and  the  Sun.  Condensation  was  already  present,  as 
the  temperature  of  the  gas  in  this  nebula  was  much  lower  than  that 
of  the  present  planets.  The  heavy  solid  particles  ranging  from 
10"^  to  10-^  cm  in  diameter  were  drawn  together  to  form  the  nucleus 
of  the  protoplanets.  This  only  represented  about  one  percent  of 
the  mass  material.  The  heavier  gases,  notably  oxygen,  nitrogen, 
the  noble  gases,  and  the  oxides  of  carbon,  nitrogen  and  sulphur, 
were  drawn  toward  the  center  and  trapped.  Many  of  the  minerals 


2 

were  carbonated,  sulphated,  and  nitrated. 

The  remaining  99  percent  of  the  mass  in  the  gaseous  envelope 
proceeded  to  evaporate  due  to  the  temperatures  generated  by  contrac- 
tion, and  due  to  solar  radiation  on  the  nearer  protoplanets .  This 
envelope  of  Hydrogen,  Helium,  ammonia  and  methane  completely  evaporated 
from  the  Terrestrial  planets.  Some  masses  broke  off  due  to  the 
centrifugal  force  of  high  rotation,  to  later  form  the  satellites  of 
Earth  and  Mars.  A  few  percent  of  the  original  envelope  remained 
with  the  Jovian  planets,  however,  and  this  accounts  for  their  re- 
ducing hydrogen  atmospheres.  Jupiter  retained  approximately  10  per- 
cent of  its  protoplanet  mass,  and  thus  is  the  largest  planet. 
Because  of  their  higher  velocities  around  the  Sun,  the  outer  planets 
lost  more  of  their  original  envelope,  thus  they  are  smaller.  It  is 
believed  that  Pluto  lost  all  of  its  gaseous  envelope  and  therefore 
closely  resembles  Earth,  except  for  its  temperature. 

Further  contraction  of  the  Terrestrial  planets  forced  percola- 
tion and  exhalation  of  gases  from  within  the  hot  core,  and  this 
accounts  for  the  Earth's  oxidizing  atmosphere.  Such  atmospheres  as 
remain  on  the  Terrestrial  planets,  then,  are  these  secondary  emissions 
of  the  heavy  gases  from  within.  Mercury  lost  most  of  its  secondary 
atmosphere  because  of  its  size  and  proximity  to  the  Sun.  Venus 
apparently  has  lost  most  of  its  water  due  to  solar  proximity,  and  Mars 
has  lost  most  of  its  atmosphere  because  of  its  small  size. 

The  evaporation  from  the  major  planets  continued  for  a  shorter 
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time  than  it  did  from  the  Terrestrial  planets.  Kuiper  estimates  10 

years  for  Earth,  and  lO^  years  for  Jupiter.  This  is  due  to  the  lower 
temperature  and  greater  distance  from  the  Sun,  and  also  to  the  larger 


mass  of  Jupiter.  The  rate  of  evat)oration  depends  orr  the  log  of 
the  mass  to  radius  ratio,  and  the  temperature.  Kuiper  has  shown 
that  Jupiter  lost  considerable  Hydrogen,  and  some  Helium  during 
the  period  of  1(K  years  when  it  could  lose  it.  At  the  present  time, 
Jupiter  is  losing  a  negligible  amount  of  gas.  We  should  note  that 
the  loss  in  H_  proceeds  as  the  square  of  the  loss  of  He,  due  to 
the  masses  of  the  two  molecules.  Therefore,  contrary  to  the 
opinion  of  Urey  (ref  2),  considerably  more  Hy  drogen  could  be 
lost  than  Helium. 

We  shall  consider  the  atmosphere  in  more  detail  in  the 
following  sections.  This  should  be  our  first  experience  with  a 
reducing,  light  gas,  "fuel"  atmosphere. 
1.2  The  Ability  to  Hold  an  Atmosphere 

According  to  the  kinetic  theory  of  gases  (e.g.,  Jeans 
(ref  3)  )»  an  atmosphere  in  approximate  isothermal  equilibrium 
has  the  following  density  profile: 

mgr/_H_\ 

/>  =  /g ." RT  r  +  H  (i-i) 

where  r  =  radius  of  planet,  H  =  height  above  datum,  and 
^  =  datum  density. 

Jeans  also  states  that  for  stability  over  astronomical 

o 
periods  (10  years),  the  RMS  molecular  velocity  of  the  gas  at 

the  exosphere  should  be  less  than  20  percent  of  the  escape 

velocity,  V^. 

V    -l/LfiL  EL.  (1.2) 

esc   (/  Rj  +  H  Sec. 

For  Jupiter,  this  works  out  to  G±.6   km/sec.  or  196,800  ft. /sec ! 


At  the  present  time,  with  temperatures  of  the  order  of  86  °K,  very  few 
molecules  at  the  escape  level  attain  this  velocity.  However,  during 
the  protoplanet  period,  with  temperatures  of  the  order  of  1000°K  in 
the  exosphere,  and  very  low  M/R  ratios,  efficient  escape  of  H?  would 
occur.  According  to  Kuiper,  the  value  of  log  M/R  for  Jupiter  at  the 
beginning  of  the  evaporative  period  was  -.39«  The  upper  limit  of 
log  M/R  permitting  efficient  escape  is:  +.18  for  m  =  2  (hydrogen), 
-.12  for  m  =  4  (helium),  and  <=*<>77   for  m  =  18  (methane,  ammonia). 
Thus  we  see  that  plenty  of  hydrogen,  some  helium,  but  no  methane 
could  have  evaporated. 

The  other  Jovian  planets  have  less  ability  to  hold  their 
envelopes,  and  it  is  estimated  that  Saturn  holds  2.5  percent  and 
Uranus  and  Neptune  about  1.25  percent.  On  this  basis  alone,  Uranus 
and  Neptune  would  be  expected  to  have  little  or  no  hydrogen,  which 
is  not  the  case.  However,  it  is  believed  that  the  extreme  cold 
temperatures  on  these  planets  may  have  slowed  diffusion  of  H2  and 
He  to  the  exosphere  during  the  evaporation  period.  Their  atmospheres 
also  show  a  lack  of  turbulence  which  is  present  on  Jupiter  and 
Saturn.  Their  relative  abundances  cannot  therefore  be  accurately 
predicted,  but  must  be  empirically  determined.  The  usual  assumption 
is  that  fractionation  factors  are  some  smooth  function  of  m,  the 
molecular  weights  of  the  constituents. 
1.3  Probable  Constituents 

The  composition  of  the  protoplanets  at  the  time  of  their 
formation  is  estimated  to  be  in  the  solar  or  cosmic  proportions. 
Accurate  spectrographic  analyses  of  the  Sun  and  stars  give  the 
relative  amounts  of  gaseous  elements  shown  in  Table  I. 


TABLE  I 

Initial  Composition  of  Gaseous 

Envelope 

Gas 

Amount 

Molecular  Weight 

Weight 

$>  Weight 

H2 

6289 

2.016 

12,679 

63.5 

He 

17^0 

4.003 

6,965 

34.9 

Ne 

5.9 

20.18 

119 

0.6 

H20 

3*75 

18.02 

68 

0.34 

NHo 

3.0 

17.03 

51 

0.26 

CH^ 

1A 

16.04 

22 

0.11 

A 

0o8 

36.3 

29 

0.15 

Totals 

8044 

19.933 

100  $ 

Average 

;  molecular  we 

ight,  m  =  2.48 

r 

=  1.46 

Urey  and  others  believe  that  the  present  remaining  envelope 
of  Jupiter  contains  hydrogen,  helium  and  methane  in  these  same 
proportions  above  the  tropopause.  To  Kuiper  it  seems  unlikely 
that  the  gases  would  remain  so  homogeneous,  and  would  escape 
uniformly,  leaving  solar  proportions  intact.  Kuiper  holds  that  the 
gases  would  be  fractioned  in  proportion  to  m.  ,  the  molecular 
weight  of  each  gas.  These  separate  opinions  lead  to  two  different 
models  of  the  atmosphere.  To  this  writer,  it  seems  more  likely 
that  Kuiper  is  correct,  and  this  theory  is  supported  by  experi- 
mental evidence  to  be  discussed  later.  However,  we  will  carry 
both  atmospheres  through  this  thesis  to  provide  a  basis  for 


comparison  of  design  parameters. 

Urey's  atmosphere,  model  "a",  is  shown  below  in  Table  II. 


TABLE  II 

Solar  Proportion  Atmosphere,  model  a 

Gas 

Amount 

Molecular  Weight 

Weight 

1*  Weight 

H 
2 

6289 

2.016 

12,679 

64.4 

He 

1740 

4.003 

6,965 

35.4 

CH. 
4 

1.4 

16.04 

22 

0.2 

Totals 

8030 

19,666 

100  i> 

Average 

molecular  weight,  m  =  2.47 

r 

■  1.47 

Kuiper's  atmosphere1;  model  b,  is  characterized  by  a  reduction 
of  Urey's  model  "a"  by  a  factor  such  that  the  ratios  of  Hydrogen  to 
Helium  to  Methane  are: 


6289n' 


H, 


1740n  :  1.4 


He  :  CH, 


(1-3) 


The  factor  n  here  is  the  reduction  due  to  evaporation,  and  as 
pointed  out  before,  Hp  evaporates  as  the  square  of  He.  To  make  the 
present  mass  come  out  right,  and  with  the  assumption  of  about  10  # 
of  the  protoplanet  envelope  remaining,  n  =  .177*  Therefore,  the 
model  b  atmosphere  is  reproduced  in  Table  III. 

We  notice  that  to  change  the  molecular  weight  by  even  this 
small  amount  requires  a  considerable  change  in  the  composition. 


TABLE  III 

Reduced 

Proportion  Atmosphere 

f  model  b 

Gas 

Amount 

Molecular  Weight 

Weight 

i  Weight 

H2 

197 

2.016 

396 

23.9 

He 

308 

4.003 

1234 

74.7 

c\ 

1.4 

16.04 

22.5 

1.4 

Totals 

506.4 

-- 

1652.5 

100. 0# 

Average 

molecular  w< 

Bight,  m  =  3.26 

1/ 

i 

=  1.56 

The  pressures  in  the  two  models  will  vary  considerably,  as  shown 
in  the  section  on  structure.  These  models  are  only  for  the  gas 
above  the  tropopause,  as  it  is  assumed  the  heavier  elements  will 
have  very  low  pressures  at  this  height.  Notice  also  that  model  b 
has  He  in  excess  of  H2  by  3  to  1,  whereas  in  the  solar  proportions 
H2  is  over  He  by  2  to  1  by  mass.  Urey  holds  that  such  a  great 
reduction  in  H2  would  be  unlikely.  Perhaps,  then,  the  actual 
atmosphere  is  somewhere  between  these  models,  for  surely  some  of 
the  gas  must  have  evaporated.  We  will  carry  both  atmospheres 
through  our  calculations  for  comparison. 
1.4  Methods  of  Measurement 

The  primary  method  used  to  determine  atmospheric  composition 
has  been  comparison  of  absorption  spectrographs  of  the  planet 
with  emission  spectrographs  of  the  Sun.  Absorption  of  a  parti- 
cular frequency  line  in  the  Sun's  spectrum  indicates  the  presence 
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of  a  particular  gas.  The  strength  of  the  absorption  indicates  the 
amount.  This  method  has  several  drawbacks,  however.  It  is  limited 
in  scope  because  the  planets  are  not  primary  emitters;  only  the  light 
from  the  Sun  can  be  re-radiated,,  Certain  gases  like  nitrogen  and 
helium  have  no  absorbtion  lines  in  the  useful  range  of  spectrography. 
The  Earth's  atmosphere  shuts  out  most  of  the  ultra  violet  wave- 
lengths. We  also  get  interference  from  CO  ,  HpO,  and  oxygen  in 
the  Earth's  atmosphere.  Most  of  these  hindrances  can  be  avoided 
with  careful  work,  but  the  most  difficult  of  all  is  to  obtain  the 
laboratory  comparison  spectra  of  the  various  gases  for  identification, 

Methane  and  ammonia  bands  were  detected  as  early  as  193^  by 
Wildt  in  the  atmosphere  of  Jupiter.  The  ammonia  band  was  much 
weaker  than  expected.  It  was  then  realized  that  the  clouds  on 
Jupiter  were  NIL  cirrus,  and  less  than  one  cm  NPT  of  free  NH~ 
existed  above  the  cloudtops.  According  to  Kuiper's  protoplanet 
theory,  the  atmosphere  should  contain  an  abundance  of  H?  and  He, 
but  no  bands  for  these  gases  were  identified.  In  1938  Herzberg 
identified  several  lines  around  7500  A°  in  the  atmospheres  of  all 
the  Jovian  planets;  and  lines  around  8270  A°  in  Uranus  and  Neptune. 
But  it  was  not  until  1952  that  a  laboratory  experiment  showed  the 
8270  A0  lines  to  be  pressure-induced  quadruple  vibration-rotation 
modes  of  the  H  molecule.  To  reproduce  this  on  Earth  required  a 
path  length  of  80  meters,  a  temperature  of  78°K  gained  by  cooling 
with  liquid  nitrogen,  and  a  pressure  of  100  atmospheres. 

This  illustrates  the  extremes  in  pressure,  path  length,  and 


cm  NPT  =  height  in  cm  of  equivalent  column  of  the  gas  on  Earth 
at  normal  pressure  and  temperature. 


temperature  that  must  be  used  to  demonstrate  on  Earth  the  same 
lines  found  in  nature's  laboratory.  It  can  be  shown  that  the 
temperature  on  Jupiter  at  high  pressures  is  too  high  to  allow 
this  excitation.  But  H2  can  be  confirmed  by  its  presence  on 
Uranus  and  Neptune,  The  lines  at  7500  A°  are  present  on 
Jupiter,  but  have  not  to  date  been  identified. 

A  rare  but  extremely  important  method  of  constitutive 
determination  was  the  occultation  of  a  star  by  Jupiter  in 
November  1952.  The  rate  of  diminution  of  light  from  the  star  as 
it  passes  behind  the  planet  is  a  function  of  the  apparent  motion, 
temperature,  surface  gravity,  and  mean  molecular  weight  of  the 
gases  in  the  atmosphere.  Since  the  first  three  of  these  parameters 
are  fairly  well  known,  the  mean  molecular  weight  can  be  found  by 
photometric  measurement  of  the  rate  of  light  diminution. 

It  was  by  this  method  that  Baum  and  Code  (ref  b)   measured 
the  m  of  Jupiter  in  1952.  The  light  curve  best  fit  an  atmosphere 
with  m  =  3«3«  This  was  important  confirmation  that  the  primary 
constituents  were  H2  and  He,  and  CH.  and  NH  as  indicated  spec- 
tographically.  It  further  indicates  an  abundance  by  weight  of 
He  over  H~  by  3«1  as  predicted  by  Kuiperl  This  is  why  this 
author  prefers  Kuiper's  model  b  atmosphere  to  Urey's  model  a. 

Other  methods  of  atmospheric  measurement  include  ultra- 
violet and  infrared  absorption  by  the  atmosphere,  polarization, 
and  radio  emission.  Although  radio  emission  in  a  band  around 
22  mc  has  been  detected  from  Jupiter,  no  conclusive  evidence  of 
composition  has  been  found.  Astronomers  believe  that  this 
radiation  is  "static",  caused  by  vast  turbulence  below  the 
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clouds  on  Jupiter.  According  to  Peek  (ref  5)  this  turbulence  gains 
its  energy  from  heat  within  the  planet,  and  from  the  fast  planetary 
rotation.  The  turbulence  accounts  for  the  meteorological  formation 
of  the  bands  on  the  planet  similar  to  the  prevailing  winds  on  Earth, 
This  will  be  discussed  more  fully  under  Structures. 
1.5  Structure  of  the  Atmosphere 

The  general  structure  of  both  models  will  conform  to  the 
following  assumptions: 

1)  A  spherically  symmetric  atmosphere 

2)  Isothermal  equilibrium  in  the  stratosphere 

3)  Dry  adiabatic  lapse  rate  in  troposphere 

*0  Moist  (NH-a)  adiabatic  lapse  rate  below  the  cloud  tops. 

From  these  assumptions,  the  thermal  models  can  be  constructed 
as  shown  in  figure  1.  The  isothermal  temperature  given  by  Peek 
(ref  5)  is  86°K.  The  equilibrium  temperature  for  solid  NH~  from 
the  phase  diagram  is  165  °K  to  l68°K,,  depending  on  pressure.  Zero 
conditions  refer  to  the  tops  of  the  clouds.  The  measurements  of 
Jupiter's  radius  have  always  been  referred  to  the  cloud  tops,  so 
we  will  use  this  as  zero  altitude. 

The  dry  adiabatic  lapse  rate  is  given  by  the  expression: 

g 

This  yields  for  Jupiter: 

1^  =  2.47)=  2.6  oK/km 

Hm  =  3.26^=  b.O  °K/km 
With  this  information,   the  depth  of  the  cloud  tops  below  the 
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tropopause  is  easily  calculated.  On  the  way  down  a  level  is 
passed  which  is  at  the  radiometric  temperature  of  Jupiter  (102°K) . 
This  is  the  level  which  the  thermocouple  sees  at  8-l^u..  It  is  also 
the  level  where  the  pressure  is  approximately  one  atmosphere*  The 
emissivity  of  the  gases  above  this  point  is  so  low  that  we  assume 
most  of  the  energy  radiated  into  space  comes  from  below  this 
level.  Between  here  and  the  cloud  tops  there  is  less  than  one 
centimeter  NPT  of  NH  ,  and  radiation  from  below  the  clouds  is 
effectively  blocked.  This  explains  the  rather  weak  ammonia 
absorption  spectrographs  from  Jupiter. 

The  density  ratio  for  the  adiabatic  portion  of  the  atmos- 
phere follows % 


=  i£_=  (ly-i)  =  (V^to 


*         T'       T 
'O  o 


(1-5) 
o 


The  density,  .^ '  ,  is  calculated  from  the  perfect  gas  law,  which 
should  apply  at  these  temperatures  and  pressures. 

p  m 

/?  = -°—  (1-6) 

r  o  r  t  v  ' 

g  o 

The  pressures  at  the  cloud  tops  are  fixed  by  the  (H2  +  He/CH.  ) 
abundance  ratio  used  for  model  a  and  b. 

p    (m  =  2.^7)  =  2h   earth  atmospheres 

P     (m  =  3- 26)  a  2  earth  atmospheres 

Here  we  point  out  again  that  model  b  is  more  probable. 
The  8270  A°  absorption  due  to  EL  quadruple  excitation  is  absent 
from  Jupiter.  But  if  a  pressure  of  2h   atm  prevailed  in  the 
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clouds  on  Jupiter.  According  to  Peek  (ref  5)  this  turbulence  gains 
its  energy  from  heat  within  the  planet,  and  from  the  fast  planetary 
rotation.  The  turbulence  accounts  for  the  meteorological  formation 
of  the  bands  on  the  planet  similar  to  the  prevailing  winds  on  Earth. 
This  will  be  discussed  more  fully  under  Structures* 
1.5  Structure  of  the  Atmosphere 

The  general  structure  of  both  models  will  conform  to  the 
following  assumptions t 

1)  A  spherically  symmetric  atmosphere 

2)  Isothermal  equilibrium  in  the  stratosphere 

3)  Dry  adiabatic  lapse  rate  in  troposphere 

*0  Moist  (NHo)  adiabatic  lapse  rate  below  the  cloud  tops. 

From  these  assumptions,  the  thermal  models  can  be  constructed 
as  shown  in  figure  1.  The  isothermal  temperature  given  by  Peek 
(ref  5)  is  86°K.  The  equilibrium  temperature  for  solid  NEL  from 
the  phase  diagram  is  165 °K  to  168°K,  depending  on  pressure.  Zero 
conditions  refer  to  the  tops  of  the  clouds.  The  measurements  of 
Jupiter's  radius  have  always  been  referred  to  the  cloud  tops,  so 
we  will  use  this  as  zero  altitude. 

The  dry  adiabatic  lapse  rate  is  given  by  the  expression: 

/?-?(*£)  tun 

© 

This  yields  for  Jupiter; 

[^  =  2.1+7)=  2.6  oK/km 

nm  =  3.26)=  4.0  °K/km 
With  this  information,  the  depth  of  the  cloud  tops  below  the 
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tropopause  is  easily  calculated.  On  the  way  down  a  level  is 
passed  which  is  at  the  radiometric  temperature  of  Jupiter  (102°K). 
This  is  the  level  which  the  thermocouple  sees  at  8-14-p,.  It  is  also 
the  level  where  the  pressure  is  approximately  one  atmosphere.  The 
emissivity  of  the  gases  above  this  point  is  so  low  that  we  assume 
most  of  the  energy  radiated  into  space  comes  from  below  this 
level.  Between  here  and  the  cloud  tops  there  is  less  than  one 
centimeter  NPT  of  NH  ,  and  radiation  from  below  the  clouds  is 
effectively  blocked.  This  explains  the  rather  weak  ammonia 
absorption  spectrographs  from  Jupiter. 

The  density  ratio  for  the  adiabatic  portion  of  the  atmos- 
phere follows '. 


C  y=l)  m   (To-Dh(ct 


j0   "  T  ;       T  K     •*' 


The  density,^  ,  is  calculated  from  the  perfect  gas  law,  which 
should  apply  at  these  temperatures  and  pressures. 

p  m 

p    = -°—  (1-6) 

/  o   R  T  K        ' 

g  o 

The  pressures  at  the  cloud  tops  are  fixed  by  the  (H?  +  He/CH.  ) 
abundance  ratio  used  for  model  a  and  b. 

P    (m  =  2.^7)  =  24  earth  atmospheres 

P     (m  =  3.26)  =a  2  earth  atmospheres 

Here  we  point  out  again  that  model  b  is  more  probable. 
The  8270  A°  absorption  due  to  H2  quadruple  excitation  is  absent 
from  Jupiter.  But  if  a  pressure  of  2h   atm  prevailed  in  the 
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visible  portion  of  the  atmosphere,  we  probably  could  see  this  band. 

Above  the  tropopause  isothermal  conditions  are  assumed  to 
prevail.  Therefore  the  density  follows  the  expression: 

f  =Ar  e  -K(H-Htr)  '  (1_7) 

Here  the  scale  height,  rp   ,  is  given  by 

a. 

gJ  m 


K  =  ^S-  (1-8) 

g 


R.Tgg 


Most  references  give  the  scale  height  of  Jupiter  as  about 
60,000  feet,  but  such  a  figure  is  based  on  a  mean  temperature  of 
the  planet  of  170°K.  Since  the  scale  height  is  to  be  used  over  the 
isothermal  portion  of  the  atmosphere,  this  author  prefers  to  use  the 
isothermal  temperature  (86°K).  This  drops  the  figure  down  consider- 
ably, and  in  fact  makes  ~  for  model  b  =  26,000  feet.  This  compares 

K 

with  the  value  for  Earth  of  23,500  feet. 

The  density  profiles  are  plotted  in  Figure  2  based  on  the  above 
assumptions.  We  should  note  that  the  density  for  model  b  is  almost 
identical  with  Earth.  Therefore,  the  dynamic  pressures  are  equal, 
and  any  winged  vehicle  which  can  fly  above  the  Earth  can  also  fly 
above  Jupiter.  Of  course  the  pressures  and  lift  required  at  a 
given  altitude  will  be  different. 
1.6  Summary 

Fairly  conclusive  evidence  has  been  offered  to  support  model  b 
over  mbdel  a.  Comparing  the  surface  of  Earth  to  the  tops  of  the 
clouds  for  Jupiter,  one  finds  the  density  profiles  nearly  identical. 
Assuming  an  entry  density  at  160  kilometers  for  Earth,  the  same 
altitude  would  prevail  for  Jupiter.  Even  with  model  a,  entry  would 
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commence  at  about  250  kilometers .  It  is  unlikely  that  any 
reasonable  shape  would  penetrate  below  the  clouds-  before  the  high 
Mach  phase  is  over.  Pressures  increase  very  rapidly  below  the 
clouds,  and  the  critical  pressure  for  Hp  is  reached  at  about 
-120  km.  From  this  point  on  we  have  solid  H  to  a  pressure  of 
200,000  atm.,  and  metalic  H  beyond.  It  seems  unlikely  that 
Earth  people  would  ever  explore  Jupiter  to  a  significant 
depth  below  the  clouds.  But  it  is  quite  likely  that  they  will 
explore  the  upper  atmosphere  for  scientific  reasons.  Jupiter  is 
the  closest  and  easiest  to  reach  of  the  remaining  protoplanets . 
The  Terrestrial  planets  all  have  secondary  atmospheres,  but  an 
accurate  determination  of  Jupiter's  atmosphere  would  give  good 
evidence  about  the  formation  of  the  solar  system. 

Such  an  entry  could  be  made  with  a  rather  unsophisticated 
probe  at  an  early  date.  The  potential  required  would  be  very 
little  more  than  the  recently  launched  Venus  probe.  Higher  g 
forces  and  heating  would  be  encountered  due  only  to  the  higher 
entry  velocities  involved.  Such  an  entry  will  be  designed  in 
the  subsequent  chapters. 
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TABLE  IV 

Summary  of  Atmospheric  Characteristics 

Property 

Jupiter  a 

Jupiter  b 

Earth 

T 

165°K 

168°K 

288°K 

p  (slugs/ft.3) 

'  0 

8.2  x  10"3 

-4 
9.02  x  10 

2.7  x  10'3 

^(atm) 

24 

2 

1 

K  (ft"1) 

2.92  x  KT5 

3.84  x  10"5 

4.26  x'10"5 

K  (km"1) 

.0958 

.1261 

5il397 

1/K  (ft) 

34,300 

26,000 

23,500 

P(oK/km) 

2.6 

4.0 

6.5 

4r  (slugs/ft3) 

1.915  x  10~3 

2.73  x  10"4 

7.07  x  1CT** 

Htr(km) 

31.5 

21.0 

11.0 

H  entry  (km) 

250 

165 

161 

r entry 

3  x  10'12 

-12 

3  x  10  ^ 

-12 
3  x  10  " 

Tiso 

86°K 

86°K 

217  °K 
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CHAPTER  2 


GENERAL  ENTRY  DYNAMICS 


2.1  Introduction 

Since  the  purpose  ~of  this  thesis  is  to  investigate  design 
parameters  for  an  unsophisticated  probe  to  Jupiter  for  atmospheric 
determination,  only  trajectories  which  lend  themselves  to  little 
or  no  guidance  will  be  considered.  This  will  limit  us  to  a  steep 
ballistic  entry  from  deep  space,  a  shallow  entry  from  a  disturbed 
reconnaissance  orbit,  or  a  decaying  orbital  entry.  The  use  of  a 
lifting  vehicle  will  be  prohibited  by  the  delicate  guidance 
required  for  orientation  and  attitude  control.  This  first  probe 
should  be  more  like  a  meteor  or  bomb,  and  may  broadcast  only  two 
things;  stagnation  point  temperature  and  specific  force.  It  may 
have  to  withstand  specific  forces  up  to  9000  Earth  "g's",  and  will 
have  to  be  dense  enough  to  make  a  deep  penetration  before  burning. 
It  will  use  a  thick  ablative  heat  shield,  and  may  rotate  slowly  to 
distribute  heat  evenly.  It  will  contain  no  gyros  or  guidance 
system,  and  its  transmitter  will  be  potted  solid  state  electronics 
and  batteries. 

In  all  these  one  pass  ballistic  entries,  the  specific  forces 
(deceleration)  will  be  far  greater  than  the  gravitational  forces 
on  the  vehicle,  and  therefore  the  entry  can  be  considered  a 
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straight  line.  A  further  assumption  is  made  of  planar  motion. 
Since  Jupiter  rotates  so  fast,  the  orientation  of  this  plane  will 
be  important.  However,  it  should  be  relatively  easy  to  make  the 
entry  near  the  equatorial  plane,  and  in  the  direction  of  rotation. 
This  will  considerably  reduce  the  entry  velocity  with  respect  to 
the  atmosphere. 

The  general  three  dimensional  equations  of  motion  for  variable 
mass;  oblate,  rotating  planet;  and  atmosphere  with  winds  are 
spelled  out  in  Duncan's  thesis  (ref  6).  These  equations  cannot  be 
solved  in  closed  form,  and  are  too  complicated  to  be  useful  in 
showing  variations  of  the  parameters  in  question.  For  a  probe  as 
straight-forward  as  this,  the  two  dimensional,  constant  mass, 
spherical  rotating  planet  without  winds  will  be  entirely  adequate. 
The  vehicle  is  essentially  a  cannonball,  and  contains  no  guidance 
to  utilize  such  nice  refinements  in  motion. 

The  ballistics  used  will  be  those  of  Chapman  (ref  7) *   and  the 
modifications  of  Duncan  (ref  6).  Duncan's  "conservation  parameter" 
will  not  be  necessary,  since  only  the  gas-dynamic  phase  is  considered. 
The  Keplerian  and  intermediate  phases  are  of  use  only  for  guidance 
in  the  range  of  10   to  10   "g's",  and  do  not  affect  vehicle  design 
for  the  final  phase. 
2.2  Types  of  Entry 

The  direct  entry  will  be  made  at  some  entry  angle ,  ^  ,  from  deep 
space  with  initial  velocity,  v  =  v   .  For  Jupiter  this  would  be 
about  200,000  ft/sec.  But  if  the  entry  were  made  in  the  direction 
of  planetary  rotation,  this  value  would  be  reduced  by  41,000 
cos  Yt   ft/sec.  This  type  entry  requires  no  rocket  thrust,  and  no 
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guidance  other  than  the  initial  aiming . 

The  shallow  entry  would  be  accomplished  by  an  impulsive 
perturbation  from  a  stable  reconnaissance  orbit  around  Jupiter. 
The  angle  would  be  very  shallow  because  large  impulsive  pertur- 
bations would  require  too  much  propellant.  A  reasonable  value 
for  ^  would  be  5Q»  The  entry  velocity  would  be  on  the  order  of 
satellite  velocity,  and  could  be  reduced  by  almost  the  entire 
rotational  velocity  of  the  planet.  Therefore  V.  =  100,000  ft. /sec. 
This  entry  would  be  gentler  than  the  first  type,  but  would  require 
a  higher  ideal  velocity  capability  for  the  launching  vehicle,  and 
much  more  accurate  guidance.  Extra  fuel  would  be  required  to 
establish  the  orbit,  and  still  more  to  launch  the  probe  from  the 
orbit . 

The  orbital  decay  entry  would  be  the  gentlest  possible  with- 
out using  a  lifting  vehicle,  but  it  is  also  the  most  difficult 
to  obtain.  A  low  orbit  within  the  sensible  atmosphere  would  have 
to  be  attained,  and  would  require  the  most  careful  guidance.  The 
time  required  for  this  entry  is  longer  than  the  others,  as  several 
passes  around  the  planet  would  be  made.  The  entry  speed  would  also 
be  about  100,000  ft/sec.  and  the  maximum  specific  force  would  be 
^0  Earth  g"s.  This  shows  that  a  lifting  vehicle  of  some  sort  will 
certainly  be  required  if  a  manned  entry  is  ever"  Made. 

At  the  present  time,  the  only  entry  available  is  the  direct 
ballistic  entry.  But  as  we  shall  show,  this  would  not  be  too 
hard  to  design  with  present  materials,  and  would  probably  give  the 
most  useful  information. 

Equations  will  be  included  for  instrumentation  of  the  probe 
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in  both  flight  path  axes  and  energy  axes.  Flight  path  axes  can  be 
used  regardless  of  the  shape  or  orientation  of  the  vehicle.  The 
assumption  has  been  made  that  for  a  steep  entry  the  specific  force 
lies  along  the  flight  path  (Drag).  Therefore,  pendulous  mounted 
instruments  will  indicate  flight  path  axes  even  in  a  spherical 
rotating  vehicle.  Energy  axes  (r,  0)  can  only  be  used,  however,  in 
conjunction  with  a  vertical  indicating  system.  The  resolution  of 
the  motion  into  radial  and  tangential  components  requires  the  angle 
)f]  •  It  is  expected  that  such  a  system  will  not  be  included  due  to 
its  complexity  and  delicate  nature. 

It  should  be  emphasized  that  this  instrumentation  is  to  deter- 
mine deceleration  forces  for  radio  relay,  and  is  not  intended  for 
guidance.  In  this  matter,  small  angle  approximation  and  neglect  of 
gravitational  terms  is  allowed. 
2.3  The  Coordinate  Frames 

Figure  3  illustrates  the  two  simple  two  dimensional  coordinate 
frames  used  in  the  analysis.  The  r,  0,  f   frame  is  a  great  circle, 
planet  centered  frame  which  rotates  with  the  vehicle.  It  can  be 
considered  an  air  mass  frame  or  energy  frame,  as  potential  energy 
and  angular  momentum  can  be  written  in  terms  of  r  and  0.  The  entry 
trajectory  is  assumed  to  lie  entirely  in  the  r,0  plane;  and  the  angle 
of  this  plane  with  the  equator,  r  ,  is  assumed  to  be  small. 

The  flight  path  coordinates  (x,z)  are  in  the  plane  of  r,0  and 
are  rotated  through  the  angle,  ¥/>  . 
2  A     The  General  Entry  Equations 
Accelerations : 

Ar  "  S  "  -T     A(4  =  ^I  +  (f  >  VI0        (2"1) 
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Components  of  Vj^: 

VI0  =  V0  +  WIJ  R  cos  If  =  V  cos  ^  +  WXJ  R  cos  f        (2-2) 

where  v  =  vel.  with  respect  to  the  air  mass, 

and  f  =  inclination  of  the  trajectory  plane  with  respect  to  the 
equatorial  plane. 
R  =  V  sin  ^  (2-3) 

It  is  useful  also  to  define  coordinates  along  and  normal  to 

the  velocity  vector  as  an  x,z  frame.  (See  Fig  3.)  The  x,z  frame 

is  rotated  from  the  r,0*  frame  by  the  angle  Zfi   „ 

Inx,z  frame; 

Ax  =  V  -  W2-^  R  cos2  f  sin  ^  (2-*0 

Az  =  V  h    -   V  °RS  6   "  ^u  R  cos2  V  cos  ^  -  2WIJ  V  cos  f 


The  forces  to  be  considered  are  the  gravitational  forces  and 
the  drag  force,, 


G  =  -  G   1 
sp  r 

=  -  G  (sin#  1  +  cos  fa   1  )  (2-5) 

sp     t>    x      b     z' 

and  drag° 

\  -  - 1  \  <2-« 

Therefore  the  two  dimensional  equations  of  motion  in  the  x,z  plane 
are; 

V  +  (G   -  W2^  R  cos2  Y  )  sin  ^  +  §  =  0   (x  force) 

(2-7) 

Vir-  V     cos  ^4-  (G       -  W2^  R  cos2n  cos<#    -  2  WTT  cos  f  =  0 
*£  R  sp  IJ  t>  IJ 

(z  force) 
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Now,  in  dimensionless  form,  normalizing  with  respect  to  RT,  V   , 

J   sat 

sat  . 

RJ 

fL+  (^  -  r^2  cos2^)  sin  £  +  ND  =  0  (2-8) 

r 

2 

v  l£  -  2-  cos  <T,  +  (i-  -  ril2  cos2!^ )  cos  ti   -   2v/lcos  V  -  0 
dT   r      »    „<£  £ 


In  r,  0  coordinates: 

dv        v  2 

^  -  -|-  +  ND  sin  <£     +  (i-  -  r/L2  cos2  ¥>)  -  2v0  Acos  f  =  0 

r  (2-9) 

S2.+  -£-£  +  ND  cos  ^  +  2vr  fL  cos  V  =  0 

These  equations  can  also  be  written  in  terms  of  energy  and 
angular  momentum,  rather  than  coordinate  forces.  Defining  zero 
potential  energy  with  respect  to  Jupiter's  gravitational  field  as 
the  energy  at  infinity,  total  energy  will  be  negative. 

j£^2ll  =  _  G  R  «.   I  (2_io) 

M        sp     2 


In  dimensionless  form: 

V2 
E  =    S(tot)   =1.1  (2-11) 

(tot)       m2  2         r 

sat 

JTJ> 

Equating  -r^.  to  power : 

^  =  f   •  v  (2-12) 

dT  I 


where  f  =  the  sum  of  all  external  forces  in  mean  surface  g's  of 
Jupiter  (dimensionless).  This  means  of  expressing  motion  is  best 
for  instrumentation,  because  accelerometers  mounted  within  the  body 
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can  then  measure  the  aerodynamic  forces,, 

-r-  =  v_  (the  output  of  the  accelerometer  in  the  1„ 
dT   Ir  r 

direction)  +  v  (the  output  of  the  accelerometer  in  the 
1,*  direction) .  (2-13) 

If  the  body  has  aerodynamic  stability,  more  useful  information 
would  come  from  body  fixed  axes,, 

S±i  =  t_  (output  of  accelerometer  parallel  to  V) 

+  v  (output  of  accelerometer  normal  to  V)        (2-1^) 

The  angular  momentum  of  the  vehicle  about  Jupiter  provides  the 
remaining  energy  expression  needed „ 
P  =  RV 

p 

non-dimensionali  zings     p  =  ~ — tt 

J     sat 
=  r   (v0  +  rAcos  *?)  (2-15) 

2£  a  r  x  f  =  r  (accelerometer  output  in  direction  lx)   (2-16) 

If  r  can  be  measured  by  radar  or  integrated  from  other  data,  the 
equations  of  motion  can  be  written  in  energy  and  angular  momentum 

form; 

2 
dE  m   d_  ,^_I   In 
dT  dT  2  ~  r 

=  (N  )(v  +  rjicos  ^  cos  f)         (energy) 

■?£:  =  -  rN_  cos  ^  (angular  momentum)  (2-17) 

dT      -U     b 

However,  to  utilize  this  system  of  instrumentation,  the  vertical 

will  have  to  be  known  to  provide  the  angle,  ^  ;  and  to  establish 
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the  r,0  coord,  system  model.  We  note  that  it  need  not  be  an  accurate 
vertical  as  would  be  necessary  for  inertial  navigation,  but  only  an 
approximation  to  the  angle,  01   .  This  could  be  obtained  from  an 
"albedo  tracker",  or  perhaps  some  other  form  of  optical  or  magnetic 
device . 

Since  refined  measurements  of  forces  are  not  needed,  it  is 
anticipated  that  the  following  data  will  be  radioed  back  to  Earth: 

1)  The  drag  load  factor,  N  ,  which  is  the  output  of  a 
pendulously  mounted  accelerometer  within  the  vehicle. 

2)  The  surface  or  internal  temperature,  or  if  the  vehicle 
does  not  rotate,  the  stagnation  point  temperature. 

2.5  General  Solutions  for  Ballistic  Entry 

These  solutions  are  all  in  non-dimensional  parameters  as 
defined  in  Appendix  I;.  They  are  based  on  Duncan's  approximations 
to  Chapman's  Z  function  analysis.  Limits  to  the  validity  of  each 
equation  are  stated  with  the  equation. 

The  basic  assumption  in  the  ballistic  entry  is  that  Oi    is 
constant.  That  is,  the.  drag  specific  forces  are  much  greater  than 
the  gravitational  forces.  This  assumption  is  valid  for  velocities 
above  25,000  ft/sec.  in  the  vicinity  of  Jupiter,  and  f or  ^  such  that: 

kX/V  >  2.5  (2-18) 

D 

For  Jupiter,  this  means  an  angle  greater  than  approximately  2.3  . 
|=  -  v  sin  ^  (2-19) 

%~t%-  -*"«**  (2-20) 

The  non-dimensional  drag  load  factor,  N^,  is  defined  as: 
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nd  •-$-**  A  Rj  ¥  <^  <2-21' 


Therefore, 

^Z  _  dv  dT_  1  Z^oR    CDS 

d<r  =  dT  do-  s  2  k  sin  ^y   M 


v  (2=22) 


Separating  variables  and  integrating,  we  have  velocity  as  a 
function  of  density  ratio,  and  atmospheric  and  vehicle  parameters s 

i         * 

Range  angle  in  the  ballistic  entry  is  found  by  the  approximations 

I  =  Xc  =  cot  #  In  £-  =  cot  <g  (h  -  h±)  (2-2^) 

It  has  been  assumed  in  the  first  chapter  that; 

<T  =  e"1*  (2=25) 

All  these  solutions  for  altitude,  density,  and  velocity  can  be 
written  in  terms  of  each  other,  or  in  terms  of  dimensionless  time, 

To 

Density  ratio  as  the  independent  variable; 

1)  h(<r)   =  -   (~)  \ncr  (2=26) 

f/°o  R       D    \  Qrtf-\ 

2)  v(<r)  =  v.    e  2k  sin£M  ~  1  (2-27) 

Altitude  as  the  independent  variable; 

1)  <T(h)  =  e~kh  (2=25) 

2)  v(h)     =  v.e\2k  sinrMAe  7  (2=28) 

x       i> 

Velocity  as  the  independent  variables 

i)  ffM  =  ^ofnc^sM  m(|-)  +  «rt  (2-29) 


ZUr 


2)       h(v)  = 


1  .        2k  sjikjTm  ,      ,v  v         -kht 
=  "  k  ln  [/oRcJs  In  (— )  +  e       * 


Dimensionless  time  as  the  independent  variable: 
1)     <TCD: 

,Ao  R   ^d^-1 

(T-  'T) 


-kv.    sin  <fcf 


e  "(2K  sin  <£ M> 


■** 


2)        h   (T): 


oo 


A  R  C_S  „  n_(j-n 


n=1         2k  sintfM' 


n  •  n  I 


•) 


-t 


ARCDS     -kh. 


^=>  a  o  R  CL.S  n 

*  k  (h.   -  h)  +    2  (-  0.      .     ")     (, 

i                 „  n  2k  sin<TM       v 

n=l  ^ 


(2-30) 
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3)      v  (T) 


/*oRCDS 
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(2-32) 


=  In 


2k^siniM  ln  (2L)  +  x 


/^oCDSR  <rt  ±n  V' 


"°    (-  ln  Mn 

n=1       n  -  nt 


+  S 


(2-33) 

The  summation  terms  in  the  functions  are  the  remaining  terms  from 
the  series  approximation  to  the  exponentials,  and  may  be  dropped  at 
the  user's  discretion. 

We  note  that  the  vehicle  parameter,  M/C  S,  turns  up  in  every 
equation,  as  well  as  the  atmospheric  parameter, /^oR. /2k,  and  the 
entry  angle  #)  .  Specific  solutions  of  the  equations  of  motion  for 
an  entry  into  Jupiter  are  presented  in  the  next  chapter.  These 
solutions  will  show  the  variation  of  the  design  parameters  for 
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specific  atmospheric  parameters  and  entry  angles, 

2,6  Solutions  for  Specific  Force  and  Stagnation  Temperature 

The  two  greatest  problems  in  planetary  entry  are  the  high 
specific  forces  and  heating  encountered,, 

The  specific  force  for  a  lifting  vehicle  at  high  Mach  is 
expressed  in  non-dimensional  form  ass 

This,  in  turn„  can  be  expressed  as  a  function  of  any  other  variable 
by  using  the  equations  in  section  2„5«  For  example,  as  a  function 
of  velocity  for  a  non-lifting  body  2 


XE  w       G£  V  2M 


a*    .   2k  sin<&M  -,_  Y_ 

1   7vv      vtJ     (2~35) 


This  expression  is  in  Earth  "g's",  and  is  valid  when  6Z  <■<  <T* 
This  function  will  have  a  maximum  when; 

3t  (vV)  "  ° 


thus 

dT   2 


2 

sin  <% 


and 


v  (L  max.)  =  o60?v„  (2-36) 

E  i 


This  result  is  independent  of  vehicle  parameters.  The 
expression  for  maximum  specific  force  depends  on  atmospheric 
and  entry  parameters  only* 

f  (max)  ^  f~-(max)  =  k(a6Q?  ^  sin  ty  (2-37) 
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The  height  at  which  this  maximum  N  occurs,  however,  is  a  function 
of  the  vehicle  design  parameter,  M/CLS. 

There  are  several  heating  factors  of  interest,  but  the  stag- 
nation point  temperature  will  be  the  hottest  to  contend  with. 
Also,  the  velocity  at  which  the  maximum  stagnation  point  temperature 
is  reached  is  independent  of  atmosphere  and  vehicle  parameters. 

Ts  =  1.392  x  104  (HF)^  (VF)1^  (v6<T)l/8  oR  (2-38) 

where  "vehicle  function"  (VF)  = =-  (2-39) 

Krad'V  RC 

and  "heating  function"  is  an  atmospheric  constant  relating  Jupiter 
to  the  Earth  as  follows: 

ratios  of: 

<j    i  w» 

(,)l/2(G)3/2(R)^(K)1/Vr)-2/3[T32L        2l/2 

~7T  J 

The  value  of  this  heating  function  for  Jupiter  is  about  64.2.  This 
will  be  discussed  more  fully  in  the  next  chapter. 

The  stagnation  temperature  as  a  function  of  velocity  is : 

T      (v)  -  1.392  x  10*  (HF)^  (VF)1^  (2*  ttjQj  v6  In  3L.)1/8 

S  EJ  fQ   R  CpS       V± 

(2-41) 

d    6 
This  is  a  maximum  when  -?=.  (v  <7")  =  0 

dv   kv  sin  <3^ 
dT      I 

and  v(T  max.)  =  .846  v.  (2-42) 

s  i 


2? 

Use  will  be  made  of  these  factors,  to  show  where  maximum 
specific  force  and  temperature  are  expected,  in  the  velocity  vs. 
altitude  plots  in  the  next  chapter. 

Little  has  been  said  of  the  decaying  orbital  trajectory. 
Such  an  entry  would  require  an  extreme  accuracy  in  guidance, 
beyond  the  scope  of  this  simple  probe.  The  esjcape  from  Earth 
and  establishment  of  a  reasonable  transfer  ellipse  (2  yr  9  mo.) 
requires  an  ideal  velocity  potential  of  about  14.3  km/sec  (ref  5)» 
To  establish  a  stable  orbit  about  Jupiter  at  all  will  require  an 
additional  ideal  velocity  of  3«7   km/sec.  This  would  be  taxing 
the  present  level  of  rocketry  to  its  limit.  A  totally  rocket 
braked  soft  landing  would  require  66   km/sec.  more,  and  is  out 
of  the  question. 

The  equations  for  establishing  a  decaying  orbit  are  out- 
lined in  Duncan  (ref  6)„  There  is  first  the  circularization 
phase  during  partial  passes  through  the  atmosphere;  then  the 
decay  phase.  These  equations  are  non-linear,  and  have  not  been 
solved  in  closed  form.  However,  machine  solutions  give  a  maximum 
ND  for  Jupiter  of  40  Earth  g's,  and  a  stagnation  point  temperature 
of  2.7  times  that  for  entry  into  Earth  (ref  7)»  Tt  appears  that 
the  unmanned  entry  is  possible,  but  the  establishment  of  the 
necessary  low  orbit  is  beyond  the  state  of  the  art  in  guidance. 
We  will  therefore  concentrate  on  the  unguided  ballistic  entry. 
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CHAPTER  3 
BALLISTIC  ENTRY  INTO  JUPITER 


3*1  Introduction 

This  chapter  deals  with  the  specific  entry  profiles  into 
Jupiter.  Solutions  to  the  equations  presented  in  the  second  chapter 
are  calculated  for  entry  into  the  two  model  atmospheres.  For  each 
atmosphere,  and  for  several  entry  angles  and  velocities,  velocity 
profiles  are  calculated.  These  show  a  variation  with  the  vehicle 
design  parameter,  M/C^S. 

Also  developed  here  are  specific  force  and  heating  profiles. 
The  heating  problem  will  be  considered  in  two  phases: 

1)  Total  heat  absorption,  Qm. 

2)  Maximum  heating  rate,  Ol^-  . 

Preliminary  calculations  have  shown  a  range  in  the  vehicle 
design  parameter  from  ^  to  100  slugs/ft.  ,  and  these  values  are 
included  in  the  figures.  Finally,  with  an  eye  toward  instrumentation, 
specific  force  and  stagnation  heating  rate  plots  against  real  time 
are  generated.  These  vary  with  design  parameter,  atmosphere  model, 
and  entry  conditions,  and  are  the  predictions  upon  which  experi- 
mental atmospheric  determination  is  suggested. 
3.2  Velocity  vs.  Altitude 

The  expression  as  developed  in  chapter  two  is  calculated  here: 
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('  °  D  y  -kh  -kh.  v 
„  ,,  v  V2k  sin^M^e  -  e  i) 
v  (h)  =  v±e      %  (3_1) 

The  plots  are  made  as  a  ratio  of  v/v.  against  altitude,  H.  There- 
fore, an  entry  from  any  initial  velocity  can  be  considered  from 
these  graphs.  The  initial  conditions  vary  with  the  model  atmos- 
phere as  follows  t 

1)  Model  a,  Figures  5s- *   b»  c 

a)  />  =  .0082  slugs/ft.3 

b)  k  =  KR 

«J  ' 

=  7615 

c)  M/CDS  =  b   to  100  slugs /ft.2 

d)  if i  =  -5°,  -30°,  -60° 

e)  h„  =  H  /R 

i    i  J 

=  .003  (\  =   500,000  ft.) 

2)  Model  b,  Figures  ^a,  b,  c 

a)  /°  =  .0009  slugs /ft.3 

b)  k  =  &mo 

c)  M/CDS  =  *J-  to  100  slugs /ft.2 

d)  ^  =  -5°,  -30°,  -60° 

e)  h±   =  .002  (H±  -  333,000  ft.) 

The  values  of  <^y  were  chosen  as  representative  of  the  direct 
entry,  and  also  the  disturbed  orbital  entry.  As  previously 
stated,  an  economical  use  of  mass  ratio  will  only  allow  the 
^  =  -5°  case  for  disturbed  orbital  entry.  These  solutions  are 
valid  to  about  .25  v. .  Beyond  this  point,  gravitational  forces 
are  significant,  and  changes  in  C^  and  other  parameters  are 
apparent  as  Mach  number  falls  below  about  ten.  However,  the 
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profiles  cover  the  range  in  which  heating  and  deceleration  effects 
are  greatest,  and  will  therefore  be  useful.  These  solutions  were 
compared  to  sample  calculations  by  machine  from  Chapman  (ref  7)* 
and  agree  closely  with  the  more  accurate  solutions  obtained  there. 

The  most  significant  finding  here  is  that  an  extremely  dense 
vehicle  will  be  required  to  reach  the  cloud  layer,  if  model  "a"  is 
accepted.  Even  at  #1   =  -60°,  M/C^S  must  be  100  to  arrive  at  the 
cloud  layer  while  still  in  the  entry  phase.  If  we  wish  to  ascertain 
the  properties  of  this  cloud  deck,  we  must  arrive  there  while  there 
is  still  a  clean  cut  response  from  either  the  heating  rate  or 
specific  force  measurements.  This  is  easily  accomplished  in  the 
model  "b"  atmosphere. 
3.3  Specific  Force 

The  specific  force  in  earth  g's  has  been  given  in  chapter  two  as: 

•p  (*t\         J  L-0   P  -J-   (rr-    x  2k  sin  oj> M  -  v  \      /0   -,c\ 
fE  (V)  =  G~    2M    V  (<7i  +  JJ   CDS   ln  v^      (2"35) 

This  force  is  assumed  to  be  all  along  the  flight  path  over  the 
portion  of  the  trajectory  where  G  is  not  considered.  In  addition, 
if  we  consider  v*  as  entry  velocity,  and  <T~.    as  negligible,  the 
maximum  specific  force  is  independent  of  vehicle  design.  Thus  the 
non-dimensional  drag  load  factor  is  given  by: 

N„  (▼)  =  J.  k  sin  42   v2  ln  Z-  (3-2) 

D      GE      *  vi 

This  load  factor  is  a  function  of  entry  angle,  entry  velocity, 
and  present  velocity.  Present  velocity,  as  a  function  of  height,  is'- 
also  a  function  of  vehicle  design.  The  maximum  load  factor,  however, 
is  a  function  only  of  the  atmosphere,  entry  angle,  and  entry  velocity, 
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The  entry  velocity  relative  to  the  air  mass  is  in  turn  a 
function  of  entry  angle  and  planetary  rotation.  In  the  case  of 
the  Terrestrial  planets,  the  rotation  is  negligible.  But  Jupiter's 
large  radius  and  fast  rotation  period  (9  hours,  50  minutes),  make 
inclusion  of  rotation  velocity  necessary.  Figure  6  shows  a-  plot 
for  both  atmospheres  of  1L.  vs  entry  angle  considering  rotation  for 
the  direct  entry.  Jupiter's  equator  lies  nearly  in  the  ecliptic 
plane,  and  we  assume  a  probe  will  remain  in  the  ecliptic  plane  for 
minimum  energy  requirements.  Therefore,  a  range  of  180°  in  entry 
angle  along  the  equator  is  covered,  considering  the  preferred 
entry  with  the  rotation,  and  the  disastrous  one  against  the  rota- 
tion. The  graph  clearly  shows  the  effect  of  the  rotational 
velocity  of  41,000  ft/sec. 

We  note  with  some  concern  that  for  the  usual  entry,  atmosphere 
b  with  the  rotation,  we  can  expect  100  earth  g's  per  degree  of 
entry  angle'.  Specially  constructed  vacuum  tubes  can  take  about 
400  g's.  Therefore  we  must  limit  the  probe  to  4°  entry  angle  if 
we  use  a  vacuum  tube  transmitter.  This  limit  is  too  restrictive 
on  the  guidance  system,  and  we  will  have  to  use  solid  state 
electronics  capable  of  about  3000  g's  for  a  30-degree  entry. 

Accelerometers  are  available  to  measure  N-  up  to  10,000  g's 
(piezo-electric  accelerometers).  The  accelerometer  system  and  its 
two  gimbal  mount  with  four  roller  bearings  will  be  designed  for 
10,000  g's.  We  note  for  interest  that  for  a  .1  pound  accelerometer, 
the  outer  gimbal  bearings  will  be  freight  car  type  roller  bearings 
with  50,000  pound  capacity  each. 

For  a  disturbed  orbital  entry,  the  maximum  N  expected  is 
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shown  in  Figure  7»  Here  the  entry  angle  is  more  easily  controlled, 
and  the  entry  velocity  is  lower „  The  expected  N_^  is  40  g's  per 
degree  entry  angle.  Also,  because  of  fuel  limitations,  a  maximum 
entry  angle  of  five  to  seven  degrees  is  expected.  The  vehicle  could 
then  use  vacuum  tube  electronics,  and  be  designed  for  250  g's.  But 
all  this  relief  comes  at  great  expense  when  it  is  realized  that  the 
vehicle  must  first  be  put  into  a  stable  orbit  about  Jupiter.  Also, 
we  see  from  recent  experience  with  the  Discoverer  satellite  series, 
that  firing  a  retro-rocket  by  remote  control  is  largely  unsuccessful. 

The  variation  in  specific  force  with  altitude  was  circulated, 
and  plotted  in  Figure  8.  This  graph  shows  how  the  vehicle  design 
parameter  controls  the  height  for  maximum  N^.  There  are  plots  for 
the  several  entry  angles  used  as  examples  for  parabolic  and  orbital 
decay  entries.  This  information  can  be  used  to  determine  the  ND 
vs  time  profile  for  the  selected  design.  It  represents  dv/dt  vs  H, 
and  therefore  H  vs  time  can  be  derived.  It  is  rather  too  much  work 
to  derive  this  for  every  possible  case,  but  it  will  be  done  for  the 
design  chosen  in  the  next  chapter.  The  specific  force  variation 
with  time  will  be  one  of  the  transmitted  outputs. 

Figure  9  is  an  illustration  of  the  time  profiles  to  be  expected. 
All  quantities  are  specified,  and  the  results  in  KL,  H,  and  V  are 
plotted.  We  see  here  that  the  entry  phase  is  very  short  compared  to 
an  entry  into  Earth.  The  deceleration  phase  takes  approximately 
three  minutes  for  Earth,  while  it  only  lasts  15  seconds  for  Jupiter. 
This  short  time  will  aid  us  in  the  structural  considerations  for  design. 
3.4  Heating 

The  expected  heat  generation  profiles  will  be  calculated  here. 
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The  problem  of  handling  the  heat  will  be  considered  with  the 
vehicle  design.  Most  sources  agree  that  for  a  ballistic  entry, 
re-radiation  of  heat  is  insignificant.  These  people  are  thinking 
in  terms  of  an  Earth  entry  of  about  three  minutes  duration.  For 
the  Jupiter  entry  of  15  seconds  duration,  there  will  be  even  less 
chance  of  a  radiative  equilibrium.  We  will,  therefore,  only 
consider  the  convective  terms  in  the  heat  flow  equations  as 
presented  by  Chapman  and  Duncan. 

This  convective  heat  will  occur  almost  entirely  in  the 
continuum  flow  regime.  According  to  the  Blasius  boundary  layer 
equation,  the  heat  transfer  through  a  laminar  boundary  layer  is 
inversely  proportional  to  the  square  root  of  the  local  Reynolds 
number  based  on  the  radius  of  curvature  of  the  nose.  Thus,  for 
bodies  in  the  same  environment,  the  heat  input  rate  per  unit  area 
and  time  will  be  inversely  proportional  to  the  square  root  of  the 
radius  curvature  of  the  missile  nose. 

In  the  ablative  heat  shield  case,  a  sharp  nosed  body  would 
soon  be  reduced  to  a  blunt  nosed  body  due  to  the  small  amount  of 
ablative  material  in  the  pointed  portion  of  the  missile.  The 
primary  reason  for  using  a  blunt  body  is  to  secure  the  high  C^ 
necessary  for  deceleration  at  high  altitude. 

The  stagnation  point  convective  heat  rate  is  given  by  the 
expression: 

(q  )  =  HgiE  (A»  (*)»  HTU  (3.3) 

c  s  Y  RC  ro   VGfl    ftr-sec 

Ca+m  is  a  function  of  the  atmosphere,  and  is  given  as  18,000  for 
Earth,  and  works  out  to  1,590,000  for  model  b  and  ^,750,000  for 


model  a.  Compared  to  Earth,  the  heating  rate  will  be  87  and  264 
times  worse,  respectively. 

In  terms  of  our  non-dimensional  parameters,  the  heating  rate: 

a  ,  Esfe  i(rf  (v)m  m —  ojh 

s   VRC  ft^-sec 

The  exponents  on  density  ratio  and  velocity  ratio  depend  on 
the  flow  conditions  in  the  boundary  layer.  This  in  turn  depends  on 
the  Reynolds  number. 

For  laminar  flow,  where  p,  is  proportional  toV  T,  then  n  ■  l/2 
and  m  =  3* 

For  turbulent  flow,  n  =  .8  and  m  =  2.2  (from  Chapman;.  Figure 
10  is  a  plot  of  Reynolds  number  at  maximum  stagnation  temperature 
velocity  (.846  v.)  for  various  M/C S  ratios.  Chapman  states  that 
we  can  have  high  hopes  of  maintaining  laminar  flow  in  the  boundary 
layer.  But  this  plot  shows  that  laminar  flow  can  hardly  be  expected 
over  much  more  than  one  foot  of  the  Jupiter  probe  vehicle.  If  we 
assume  the  boundary  layer  will  become  turbulent  at  6  or  7  million, 
a  body  of  M/C-pS  =  10  at  an  entry  angle  of  30  degrees  will  have 
laminar  flow  for  about  one  foot  from  its  leading  edge.  However, 
if  the  body  is  blunt  enough  (e.g.,  like  a  Mercury  Capsule)  most  of 
the  heating  will  occur  over  the  laminar  flow  regime. 

These  extremely  blunt  bodies  will  suffer  more  from  radiation 
heating  from  the  hot  gas  to  the  vehicle,  but  the  time  of  flight  is 
so  short  we  believe  radiation  will  have  little  effect  on  the  total 
heating  problem. 

The  maximum  laminar  heating  rate  is  given  by  the  expression: 

(Q  )  (max)  =  ^  (.605)  (<T)l/Z   (v.  )3  (3-5) 

c  s      Vrc  1 
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This  is  a  function  of  M/C  S  and  entry  conditions.  The  density- 
ratio  used  is  that  for  which  v  =  .8^6  v  ,  and  is  calculated  from 
the  graphs  in  Figures  k-  and  5»  Figure  11  is  a  plot  of  the  maximum 
laminar  stagnation  point  heating  rate  expected  vs.  the  design 
parameter  M/C~S.  Also  shown  on  this  graph  for  comparison  are  the 
maximum  allowable  heating  rates  for  Teflon  and  graphite  ablative 
heat  shields  on  a  spherical  vehicle  (ref  8,9).  This  graph 
indicates  at  once  that  we  will  not  be  able  to  tolerate  steep 
ballistic  entries  with  our  present  ablative  protection.  This  is 
one  of  the  important  constraints  to  be  used  in  the  next  chapter. 

Figure  12  is  an  example  of  a  heat  rate  vs.  time  plot  for  a 
particular  entry  configuration.  The  integral  of  plots  such  as 
these  will  give  the  total  heat  absorbed,  and  will  determine  the 
weight  of  ablative  material  to  be  used  in  the  vehicle  design. 
3.5  Summary 

We  have  now  developed  particular  constraints  for  the  Jupiter 
probe.  The  structural  limitations  under  high  specific  force*  and 
the  heat  dissipation  problem  will  determine  the  upper  limit  of 
severity  of  the  entry.  Fortunately,  human  factors  are  not  involved, 
The  increase  in  guidance  requirements,  and  mass  ratio  of  the 
booster  will  determine  the  practical  floor  on  the  gentleness  of 
the  entry.  In  the  next  chapter,  we  shall  see  whether  these 
constraints  leave  a  channel  of  possibility  for  the  Jupiter  probe 
using  present  day  hardware. 

The  heating  rates  developed  here  are  for  a  laminar  boundary 
layer.  The  heating  rates  in  a  turbulent  boundary  layer,  due  to 
better  convection,  would  be  prohibitive.  But  we  must  not  fool 
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ourselves  into  thinking  the  boundary  layer  will  remain  laminar  very 
far  aft  of  the  nose.  This  requirement  on  boundary  layer  will  present 
a  further  constraint. 

Stagnation  point  temperature,  per  se,  has  little  relation  to 
the  entry  profile  when  an  ablative  coating  is  used.  The  shield 
quickly  heats  to  the  vaporization  temperature,  and  remains  constant 
while  the  heat  is  dissipated.  This  boiling  off  of  the  gaseous 
material  further  shields  the  vehicle  itself  from  convective  heat 
input.  The  vehicle  skin  will  remain  below  the  critical  2000°F 
temperature  (for  steel)  as  long  as  the  coating  lasts.  This  is  due 
primarily  to  the  poor  heat  conduction  of  the  ablative  coating. 
These  considerations  will  influence  the  instrumentation  of  the  probe. 

The  problem  of  communication  of  the  data  back  to  Earth  is  not 
considered  in  this  thesis.  However,  it  should  be  recognized  that  an 
ionized  shock  layer  may  interfere  with  electromagnetic  transmission 
at  certain  frequencies.  The  detail  design  of  the  antenna,  and  the 
selection  of  transmitting  frequency,  will  have  to  include  considera- 
tion of  the  degree  of  ionization  in  the  atmosphere  around  the  vehicle. 
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CHAPTER  4 

CONSTRAINTS  ON  VEHICLE  DESIGN  AND  FEASIBILITY 
OF  ATMOSPHERIC  DETERMINATION 


4.1  Introduction 

At  this  point,  we  must  examine  the  constraints  imposed  on  the 
vehicle  design  by  the  forces  and  heating  rates  calculated  in  the 
previous  chapters.  We  will  then  estimate  the  constraints  due  to 
the  present  limitations  in  guidance.  The  vehicle  configuration 
will  be  specified.  Within  the  area  of  possibility  for  a  success- 
ful entry,  we  will  discuss  the  feasibility  of  atmospheric  model 
determination.  The  conclusions  and  recommendations  will  follow. 

4.2  Heating  Rate  Constraints 

We  have  discovered  that  the  only  feasible  means  of  heat 
dissipation  will  be  an  ablative  shield.  In  view  of  this,  and  the 
short  time  duration  of  the  heat  input  (2  to  10  seconds),  the  heating 
of  the  vehicle  core  should  be  acceptable.  We  can  carry  enough 
heat  shield  without  a  great  payload  penalty  to  handle  even  this 
large  amount  of  heat.  The  vaporization  of  the  material  will  not 
only  carry  away  heat  to  the  surrounding  air,  but  will  provide  an 
insulating  region  of  vaporized  material  to  further  shield  the 
vehicle  from  heat  convection.  Despite  the  large  heat  inputs,  the 
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materials  used  are  poor  conductors,  and  the  interior  should  remain 
below  its  critical  strength  temperature  (2000°F  for  steel) . 

The  material  most  generally  used  for  this  type  shield  is  Teflon 
(a  Dupont  trademark  for  polyethof lourethylene) .  It  is  a  plastic 
with  strength  at  high  temperatures,  and  is  easily  molded.  Its  heat 
capacity  to  vaporize  is  about  10,000  BTU/lb.,  and  the  maximum  heat 
rate  input  is  about  50,000  BTU  /ft.2sec.  (ref.  8).  This  limit  keeps 
the  material  from  boiling  off  so  fast  that  it  explodes. 

The  best  material  for  an  ablative  heat  shield  from  the  point  of 
view  of  heat  capacity  is  graphite.  Graphite  can  dissipate  about 
30,000  BTU/lbs.,  and  at  a  rate  of  150,000  BTU/ft.2sec.  But  there 
are  many  structural  problems  associated  with  the  use  of  graphite 
that  have  not  yet  been  solved.  It  is  brittle,  and  hard  to  cast 
unless  it  is  mixed  with  other  more  plastic  compounds.  But  such  a 
mixtvre  lowers  the  heat  capacity.  Perhaps  a  mixture  of  Teflon  and 
graphite  may  be  successfully  used. 

The  heat  rate  input  graph  (Figure  11)  shows  the  limits  imposed 
by  these  materials  on  the  type  of  possible  entry.  The  Teflon 
coated  sphere  of  reasonable  density  (from  1  to  h  with  respect  to 
water)  can  only  be  used  for  the  deflected  orbital  entry.  If  we 
allow  a  solution  to  the  graphite  structural  problem,  we  can  still 
only  make  a  20  degree  direct  entry.  These  are  the  rather  stringent 
constraints  imposed  by  the  maximum  allowable  heat  rates. 
^-.3  Configuration 

Another  glance  at  Figure  11  shows  us  that  the  lighter  vehicles 
will  allow  a  steeper  trajectory  than  the  heavy  ones.  Most  of  the 
missiles  used  today  have  a  payload  specific  weight  very  close  to 
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one  (with  respect  to  water) .  Figure  13  shows  the  variation  in 
M/CpS  with  radius  of  curvature  for  a  sphere.  Figure  1*4-  shows  this 
information  for  a  typical  final  booster  stage,  and  Figure  15*  for 
a  project  Mercury  shaped  capsule.  A  C  of  1.0  is  used  for  these 
shapes  a"fc  hypersonic  velocities.  The  Pioneer  V  probe  has  shown  a 
deep  space  launching  capability  of  about  100  pounds  payload.  The 
Saturn  booster  is  reported  to  have  a  15,000  pound  payload.  These 
lines  are  also  indicated  on  the  configuration  graphs. 

It  is  apparent  that  if  we  limit  ourselves  to  100  pounds,  we 
can  only  attain  M/C  S  of  about  4  slugs /ft.  .  This  would  be  to  our 
advantage  for  heat  protection,  but  would  not  allow  penetration  to 
the  vicinity  of  the  cloud  level  during  the  entry  phase.  Using  the 
Saturn  booster,  we  can  attain  M/C_S  on  the  order  of  20  slugs/ft. 
This  would  give  us  the  desired  profile  in  atmosphere  model  b.  We 
would  have  to  go  to  still  greater  weights  to  reach  the  clouds  for 
atmosphere  model  a. 

We  note  that  the  project  Mercury  qapsule  has  the  advantage  of 
a  larger  radius  of  curvature  for  its  diameter  than  the  other 
typical  shapes  considered.  This  shape  has  been  demonstrated  to 
have  aerodynamic  stability  at  high  Mach  number,  but  a  still  flatter 
nose  would  probably  lead  to  instability.  In  Figure  15  we  have  run 
the  calculations  to  a  density  of  7»  This  would  be  a  solid  steel 
vehicle  with  a  graphite  heat  shield,  and  would  approximate  the 
limit  of  density  in  vehicle  design. 

In  order  to  show  variation  of  design  to  best  advantage  in 
terms  of  mission,  entry  angle,  and  M/C~S  we  must  assume  some 
basic  limits.  We  assume  that  the  Saturn  booster  is  available, 
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and  that  a  graphite  heat  shield  can  be  used.  We  therefore  let  the 
weight  equal  15,000  pounds.  We  allow  average  vehicle  density  to  run 
from  7  to  0.3  with  respect  to  water.  We  limit  the  vehicle  size  to 
a  10-foot  diameter. 

Using  these  criteria,  and  satisfying  the  maximum  heat  rate 
constraint,  we  develop  the  plot  in  Figure  16.  The  variation  of  M/C^S 
with  entry  angle  as  a  function  of  vehicle  average  density  for  the 
three  shapes  chosen  is  shown.  The  mission  requirement  is  arbitrarily 
set  as  reaching  the  clouds  with  v  =  .25  v. .  From  Figure  16,  several 
facts  are  evident. 

1)  With  the  very  low  M/C  S  compatible  with  the  present  100 
pound  launch  capability,  the  vehicle  cannot  hope  to  reach  low 
altitudes  during  the  entry  phase.  We  must  therefore  wait  for  the 
Saturn  booster. 

2)  The  corlifcal  shaped  vehicle  has  an  advantage  due  to  its 
larger  radius  of  curvature  for  a  given  vehicle  diameter.  This 
vehicle  can  come  fairly  close  to  meeting  the  mission  requirement. 

3)  It  would  tax  the  imagination  to  conceive  of  a  vehicle 
large  enough  to  penetrate  atmosphere  model  a  to  the  cloud  level. 
M/CjJS  would  have  to  be  100  slugs /ft.  ,  and  with  a  solid  steel  vehicle, 
the  diameter  would  be  20  feet,  and  the  weight  would  be  127,000  pounds. 
In  addition,  the  ablative  material  would  have  to  absorb  ^50,000 
BTU/ft.2  sec. 

k)     There  will  be  a  lower  limit  on  the  vehicle  average  density 
dictated  by  the  specific  force  requirements.  A  detailed  structural 
analysis  is  beyond  the  scope  of  this  thesis.  The  specific  force 
goes  up  practically  linearly  with  entry  angle  if  we  enter  in  the 


direction  of  planetary  rotation,  and  figures  to  about  100  Earth 
g's  per  degree.  We  suggest  that  a  reasonable  structural  limit  for 
a  design  load  factor  of  3000  would  produce  an  average  vehicle 
density  of  two. 

5)  It  is  evident  from  the  figure  that  high  M/C_S  alone  does 
not  guarantee  a  deep  penetration.  The  heating  rate  constraint 
limits  such  dense  vehicles  to  a  low  entry  angle  for  parabolic 
entry  velocity. 

6)  The  deflected  orbital  entry  is  virtually  free  from  heating 
constraints.  Any  of  these  vehicles  with  an  M/C-.S  of  100'  can  enter 
at  more  than  30°  without  overheating.  But  the  deflected  orbital 
entry  is  not  that  easily  accomplished.  We  are  limited  by  fuel 
carried  in  the  vehicle  to  a  deflection  of  5  or  10  degrees.  From 
Figure  16,  we  must  have  a  high  M/CLS.  At  5°  deflection  we  must 
have  100  slugs /ft.   Only  the  rocket  shaped  vehicle  can  reach 

that  high  an  M/C~S  without  exceeding  the  15.000  pound  weight  limit. 

The  configuration  chosen  as  most  likely  to  succeed  is  the 
Mercury  capsule  shape  with  an  average  density  of  two.  The  entry 
angle  should  be  near  30°.  The  M/C~S  is  11,  the  nose  radius  and 
vehicle  diameter  is  8.5  feet,  and  the  weight  is  15»000  pounds. 
A  graphite  heat  shield  is  used,  and  the  maximum  heating  rate  will 
be  150,000  BTU/ft.2sec.  The  significant  heating  period  is  2  A 

seconds,  and  by  graphical  integration  the  total  heat  absorbed  is 

2 
180,000  BTU/ft.  of  heat  shield  (Figure  12).  This  specifies  the 

1 
0 
shield  weight  as  6  pounds /ft.   (Ref .  9).  and  makes  total  ablative 

material  weight  I36O  pounds,  or  about  9$  of  the  vehicle  weight. 

Other  materials  or  cooling  systems  would  require  much  greater 
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weight,  and  this  9$  seems  to  be  a  practical  low  for  the  Jupiter  entry 
vehicle.  The  interior  will  stay  relatively  cool  as  long  as  there 
is  some  of  the  ablative  shield  remaining.  If  we  make  the  shield 
lighter,  the  structure  would  melt  or  lose  strength  as  soon  as  the 
ablative  material  evaporates.  The  majority  of  the  deceleration 
loads  occur  after  the  heating  phase,  and  we  must  not  allow  the  vehicle 
to  become  weakened  early  in  the  trajectory. 

The  heating  protection  described  here  is  for  the  model  b 
atmosphere.  From  Figure  11  it  is  evident  that  only  a  5°  entry 
could  possibly  be  tolerated  in  atmosphere  a.  Such  an  entry  would  not 
meet  mission  requirements,  or  guidance  constraints.  It  is  only 
practical,  then,  to ^design  the  vehicle  for  model  b,  and  allow  its 
failure  in  model  a  to  be  an  indication  of  the  atmosphere. 
4.4  Guidance  Constraints 

It  is  difficult  to  say  just  how  good  our  present  guidance 
systems  are.  The  Pioneer  V  probe  to  the  orbit  of  Venus  was  judged 
to  have  performed  perfectly.  Yet  it  is  short  of  its  goal  by  several 
thousand  miles  due  to  a  booster  burning  time  uncertainty  of  less  than 
one  second.  It  seems  clear,  then,  for  a  precision  injection  at  any 
controlled  angle  into  Jupiter,  we  must  have  mid-course  and  terminal 
guidance.  Such  a  system,  furthermore,  must  operate  automatically 
without  signals  from  Earth.  It  must  take  its  measurements  in  the 
vicinity  of  Jupiter  from  the  fields  of  Jupiter  (optical,  magnetic, 
or  electric).  It  would  probably  carry  a  model  inertial  reference 
frame  (stable  platform  with  gyros),  which  would  be  corrected  for 
drift  periodically  by  star  trackers  during  its  2  year  9  month  journey. 

This  system  does  not  exist  today,  as  far  as  the  unclassified 
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literature  is  concerned „  If  it  did,  it  would  have  been  used  in  the 
Pioneer  V  probe  to  Venus.  Such  developments  may  reasonably  be 
expected  in  the  next  ten  years,  however.  The  consensus  at  this 
Institute  is  that  with  our  present  launch  guidance  we  would  be 
lucky  to  hit  Jupiter  at  all.  The  Keplerian  trajectory  and  ballistic 
entry  make  no  provision  for  corrections  to  launch  errors,  and 
certainly  we  could  not  expect  a  precise  injection  angle,  much  less 
a  stable  orbit  about  Jupiter. 

It  might  be  interesting  to  consider  the  statistical  chances 
of  attaining  a  30°  or  less  injection  in  the  direction  of  planetary 
rotation,  assuming  we  can  attain  a  parabolic  contact  with  the 
planet.  This  would  be  an  academic  study  at  best,  for  we  must  not 
expect  the  United  States  to  expend  this  awful  amount  of  money 
without  a  good  probability  for  success.  It  seems  we  are  con- 
strained to  wait  for  a  system  including  terminal  guidance. 
4.5  Feasibility  of  Atmospheric  Determination 

For  this  discussion,  we  shall  assume  we  can  make  the  30° 
entry  with  the  vehicle  selected  under  section  4.3.  What  informa- 
tiorJshould  we  instrument  and  radio  back  to  Earth? 

An  examination  of  Figure  6  indicates  that  there  is  considerable 
difference  in  the  maximum  specific  force  expected  in  the  two 
atmospheres.  Model  a  will  present  an  N^  (max)  2l\7% lower  than 
model  b.  This  difference  is  quite  sensitive  to  entry  angle,  how- 
ever, and  an  error  of  10°  in  entry  angle  will  also  equal  24$  of 
N-p..  We  must  therefore  include  instrumentation  to  measure  the 
entry  angle.  Even  if  this  angle  cannot  be  controlled,  the  model 
atmosphere  can  be  ascertained  if  the  angle  and  its  degree  of 
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certainty  are  known.  The  maximum  specific  force  is  a  definite 
quantity  not  dependent  on  the  vehicle  configuration,  and  therefore 
not  subject  to  uncertainties  such  as  CL.,  or  effective  drag  area 
(which  may  be  considerably  larger  than  the  cross  sectional  area  of 
a  vehicle  in  an  ionized  gas).  This  force  is  a  function  solely 
determined  by  the  entry  parameters  (which  we  have  assumed  ar© 
successfully  controlled),  and  the  atmospheric  parameter,  k,  to  be 
determined . 

The  skin  temperature  or  its  derivative  appear  to  hold  atmospheric 
information  for  us  also.  Figure  10  shows  that  the  Reynolds  numbers 
for  the  two  models  are  nearly  identical.  However,  the  heating 
profiles  are  quite  different.  The  higher  cloud  level  pressure  and 
density  of  model  a  account  for  the  higher  heating  rates  (Figure  11) . 
Unfortunately,  a  temperature  would  be  difficult  to  instrument  in  a 
vehicle  incorporating  an  ablative  heat  shield.  The  heating  is  wholly 
transient,  and  in  the  short  effective  heating  period,  no  equili- 
brium temperatures  are  reached.  Further,  the  temperature  or  heating 
rate  are  quite  sensitive  to  uncertainties  such  as  C_.,  S,  Prandtl 
number  and  viscosity.  To  use  heating  as  a  determinant,  these 
instrumentation  difficulties  and  uncertainties  would  have  to  be 
overcome . 

A  very  useful  determinant  would  be  altitude,  and  its  measure- 
ment by  radar  would  also  provide  the  entry  angle  K   ,  needed  for  N^ 
determination.  The  height  at  which  maximum  Np  occurs  varies 
considerably  with  the  atmospheric  parameter  (Figure  8) .  The  clouds 
of  Jupiter  are  expected  to  be  ammonia  cirrus  (solid  crystals),  and 
would  give  a  good  radar  echo  to  the  proper  frequency  emission.  The 
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radar  could  be  made  to  track  minimum  range*  and  the  dish  angle 
would  then  be  the  flight  path  angle. 

There  would  be  difficulties  with  the  radar,  of  course.  The 
tracking  mechanism  would  have  to  be  perfectly  balanced  and  very 
rugged  to  endure  3000  g's.  The  electronics  would  have  to  be  all 
solid  state.  There  may  be  interference,  as  with  the  radio  trans- 
mission, from  the  ionized  shock  layer.  But  these  problems  are  not 
insurmountable,  and  should  require  no  more  development  time  than 
the  guidance  system  required. 

Even  without  the  radar,  specific  force  measurements  would  be 
indicative  of  the  atmosphere  to  the  degree  of  uncertainty  of  the 
entry  angle. 
^-.6  Conclusions  and  Recommendations 

The  major  problem  appears  to  be  control  of  the  entry  angle. 
Atmospheric  determination  can  be  made  if  the  entry  is  successful, 
and  if  the  entry  angle  is  known  within  about  five  degrees.  To 
make  a  successful  entry,  the  entry  angle  must  be  controlled  to 
30°  or  less. 

The  mission  is  stymied  at  present  by  the  lack  of  a  mid- 
course  and  terminal  guidance  system  that  would  provide  this 
accuracy  in  entry  angle.  Such  a  system  is  probably  under  develop- 
ment at  this  time.  It  need  not  be  stressed  for  the  high  load 
factors  encountered  during  the  entry  phase.  It  would  be  sufficient 
to  control  the  entry  angle  prior  to  the  onset  of  the  specific  forces. 

If  this  controlled  entry  angle  were  reliable,  the  only 
instrumentation  needed  in  the  entry  vehicle  would  be  a  pendulously 
mounted  piezo-electric  accelerometer.  A  potted  electronics 


transmitter  could  be  built  to  broadcast  the  load  factor  thus 
measured.  If  the  controlled  entry  were  not  reliable  to  within 
five  degrees,  a  device  would  be  needed  in  the  probe  to  measure  the 
actual  angle.  If  radar  is  developed  for  this  measurement,  it  can 
also  supply  the  vehicle  height  above  the  clouds.  This  height  can 
be  used  as  a  cross  check  in  atmospheric  determination. 

The  vehicle  itself  will  have  to  be  large.  In  the  Mercury 
capsule  configuration  the  vehicle  could  reach  nearly  to  the  cloud 
level  at  a  weight  of  15,000  pounds.  This  will  require  the  Saturn 
booster  for  launch  and  a  graphite  heat  shield  —  both  under  develop- 
ment at  this  time.  The  present  deep  space  launch  capability  of  100 
pounds  will  not  be  adequate  for  the  Jupiter  probe.  A  Teflon  heat 
shield  does  not  have  enough  heat  rate  capacity  for  the  direct  entry, 
but  would  do  for  the  deflected  orbital  entry. 

The  decaying  orbital  and  deflected  orbital  entries  are  still 
farther  into  the  future  with  respect  to  guidance. 

The  primary  purpose  of  the  mission  is  to  determine  the  molecular 
weight  of  Jupiter's  atmosphere.  It  has  been  shown  that  this  can 
be  done  with  a  large,  but  simple,  entry  vehicle.  But  the  entry 
vehicle  must  be  delivered  to  the  entry  point  by  a  much  more  sophis- 
ticated navigational  vehicle.  This  delivery  is  beyond  the  present 
state  of  the  art  of  guidance. 

We  suggest  for  further  study  the  following  points: 

1)  The  detailed  structural  design  of  the  vehicle  to  withstand 
the  high  load  factors  encountered; 

2)  The  degree  of  ionization  in  the  shock  layer,  and  its  effect 
on  electro-magnetic  transmission; 
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3)  The  chemical  effect  of  the  heat  and  the  reducing  atmos- 
phere on  the  ablating  material; 

k)     The  probability  of  achieving  a  30°  or  less  entry  angle 
with  present  day  launch  guidance  only. 

We  believe  the  guidance,  booster,  and  graphite  heat  shield 
problems  are  currently  under  consideration. 
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A .  Subscripts 

i 
I 
J 
E 
b 
s 

B .  Symbols 

A 


WD 

°V 
D 

e 

E 

F 
t 


APPENDIX  I 
NOTATION  AND  SYMBOLS 

Components  along  the  ly,  ,  1^  ,  1^  axes. 

Initial  value  of  the  quantity. 

Measured  with  respect  to  inertial  coordinates. 

The  quantity  for  Jupiter. 

The  quantity  for  Earth. 

Refers  to  ballistic  entry. 

Refers  to  stagnation  area. 

Acceleration  vector  of  entry  vehicle,  ft/sec  . 


C  .   Dimensional  constant  for  the  planetary  atmosphere. 

3/2    -1 
Drag  coefficient.  (BTU  ft    seC  > 


Heat  capacity  at  constant  pressure. 
Heat  capacity  at  constant  volume. 
Drag  force  (lb.). 
Energy  (ft-lb) . 
Dimensionless  energy  =  — " 


MV* 


sat 


Specific  force  (ft /sec.  ) 


Dimensionless  specific  force  in  surface  g's  of 


the  planet. 


f  = 


G 


G  Gravitational  field  intensity  (ft/sec  ) 

—  2 

g  Planet's  gravity  field  (ft/sec  ) 

H  Altitude  of  vehicle  above  planetary  reference  level  (ft) 
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h     Dimensionless  altitude  =  -5— 

A 

K     Exponential  decay  parameter  of  planetary  atmosphere 

(ft'1). 
K .  .   Vehicle  surface  radiation  emissivity. 

k     Dimensionless  decay  parameter  of  planetary 

atmosphere  =  KR 
M     Mass  of  vehicle  (slugs) 
M.     Mass  of  planet  (Jupiter,  in  this  case) 
Nn     Dimensionless  drag  load  factor  -  r^— 

JJ  Ma 

P     Angular  momentum  per  unit  mass  of  the  vehicle 

p 

p      Dimensionless  angular  momentum  =  — ry= — r-rr- 

R    G 

P^     Prandtl  number 
r 

Q     Total  convective  heat  absorbed  per  unit  area. 

ft 

RC     Radius  of  curvature  of  vehicle  nose  (ft). 
R     Universal  gas  constant  =  8. 31  x  10'  ergs/°K 

o 

R     Radius  vector  from  planet  center  to  vehicle 

T  Dimensionless  radius  vector  =  - — 

RJ 
Re     Reynolds  number  =/?V  lenKth 

S     Reference  area  of  entry  vehicle  used  in  drag 

computations  (ft  ). 
T      Temperature  (°K). 
t      Real  time  (sec). 
V     Velocity  vector  of  a  vehicle  with  respect  to 

coordinates  rotating  with  a  planet. 
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Dimensionless  velocity  = 


sat 


V  i 

sat  Circular  orbital  velocity  at  planet  surface  =  y    GR 

X    Distance  flown  measured  at  surface  of  planet 
Z     Chapman's  transformation  variable  =  v  cos  *b  sin  *b 

k 


C.  Greek  Symbols 


y  Ratio  of  specific  heats  C D/C 

y^  Ballistic  entry  angle 

Jf  Dry  adiabatic  lapse  rate  (°K/km) 

\x  Coefficient  of  viscosity  of  the  planet's  atmosphere 

(slugs /ft- sec) 

p  Free  stream  atmospheric  density  (slugs /ft-') 

jO  Reference  level  atmospheric  density 

<T  Density  ratio   /° l/> 

S(  )  Summation  of  quantities  in  parenthesis 
n=l 

T*    Dimensionless  time  -    — g —  t 

0    Angle  measured  in  plane  of  trajectory  from  a 
reference  point  in  the  direction  of  motion 

y    Angle  of  inclination  of  the  trajectory  plane  with 
the  planet's  equatorial  plane. 

G     Dimensionless  angular  velocity  of  the  planet  about 
its  polar  axis 


ij/g 


Q  =  ^jt/g     for  ^piter 
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1±   initial  range  angle 


FIGURE  3  ENTRY  COORDINATE  SYSTEMS 
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FIGURE  11  MAXIMUM  STAGNATION  AREA  HEATI 


vs  M/CpS 
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